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ABSTRACT 

When approximating a space curve, it is natural to consider whether the knot type 
of the original curve is preserved in the approximant. This preservation is of strong 
contemporary interest in computer graphics and visualization. We establish a criterion 
to preserve knot type under approximation that relies upon pointwise convergence and 
convergence in total curvature. 
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1. Introduction 

Curve approximation has a rich history, where the Weierstrass Approximation The- 
orem is a classical, seminal result |23| . Curve approximation algorithms typically 
do not include any guarantees about retaining topological characteristics, such as 
ambient isotopic equivalence. One may easily obtain a sequence of non-trivial knots 
converging pointwise to a circle, with the knotted portions of the sequence becom- 
ing smaller and smaller. These non-trivial knots will never be ambient isotopic to 
the circle. However, ambient isotopic equivalence is a fundamental concern in knot 
theory. Moreover, it is a theoretical foundation for curve approximation algorithms 
in computer graphics and visualization. 

So a natural question is what criterion will guarantee ambient isotopic equiva- 
lence for curve approximation? The answer is that, besides pointwise convergence, 
an additional hypothesis of convergence in total curvature will be sufficient, as we 
shall prove. An example is shown by Figure [l] 



Figure 1(a) shows a knotted curve (yellow) which is a trefoil, where this curve is a 
spline initially defined by an unknotted PL curve (purple), called a control polygon. 
This PL curve is often treated as the initial approximation of the spline curve. A 
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(a) Unknot vs. Knot (b) Knot vs. Knot 

Fig. 1. Ambient isotopic approximation 

standard algorithm, called subdivision [S', is used to generate new PL curves that 
more closely approximate the spline curve. Figure |l(b) shows an ambient isotopic 



approximation generated by subdivision, as this PL approximation is a trefoil. 

There are three main theorems presented. All have a hypothesis of a sequence 
of curves converging to another smooth curve C. In Theorem |4.6[ the elements of 
the sequence are PL inscribed curves. In Theorem |5.3| and |7.8[ the class of curves is 
generalized to any piecewise curves, with the first being a technical result about 
a lower bound for the total curvature of elements in some tail of the sequence. 
These first two results are used to provide the main result Theorem |7.8[ showing 
that pointwise convergence and convergence in total curvature over this richer class 
of piecewise curves produce a tail of elements that are ambient isotopic to C. 



2. Related Work 

The Isotopic Convergence Theorem presented here is motivated by the question 
about topological integrity of geometric models in computer graphics and visual- 
ization. But it is a general and pure theoretical result, dealing with the fundamental 
equivalence relation in knot theory, which may be applied, but extends beyond the 
limit of any specific applications. 

The preservation of topology in computer graphics and visualization has previ- 
ously been articulated in two primary applications 9J: 

(1) preservation of isotopic equivalence by approximations; and 

(2) preservation of isotopic equivalence during dynamic changes, such as protein 
unfolding. 

The publications |1|2)15|18] are among the first that provided algorithms to en- 
sure an ambient isotopic approximation. The paper [Mi provided existence criteria 
for a PL approximation of a rational spline curve, but did not include any specific 
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algorithms. 

Recent progress was made for the class of Bezier curves, by providing stopping 
criteria for subdivision algorithms to ensure ambient isotopic equivalence for Bezier 
curves of any degree n [10], extending the previous work of jl8], that had been re- 
stricted to degree less than 4. This extension is based on theorems and sophisticated 
techniques on knot structures. 

This work here extends to a much broader class of curves, piecewise curves, 
where there is no restriction on approximation algorithms. Because of its general- 
ity, this pure mathematical result is potentially applicable to both theoretical and 
practical areas. 

There exist results in the literature showing ambient isotopy from a different 
point of view |4I24| . Precisely, there is an upper bound on distance and an upper 
bound on angles between corresponding points for two curves. If the corresponding 
distances and angles are within the upper bounds, then they are ambient isotopic. 

Milnor [T^ defined the total curvature for a curve using inscribed PL curves. 
The extension of the definition to piecewise curves can be trivially done. Con- 
sequently, Fenchel's Theorem can be applied to piecewise curves, as we need 
here. 

Milnor [IB] also proved the ambient isotopy between a given curve and the 
inscribed curves. That is a similar version of Theorem |4.6| presented here. That 
result was recently generalized to finite total curvature knots The application 
to graphs was also established recently [7] . Our proof here indicates an upper bound 
on distance and an upper bound on total curvature for ensuring the isotopy, which 
leads to the formulation of algorithms. 



3. Preliminaries 

Use C to denote a compact, regular, C^, simple, parametric, space curve. 
Let {Ci}^ denote a sequence of piecevifise C^, parametric curves. Suppose 
all curves are parametrized on [0,1], that is, C — C{t) and Ci — Ci{t) for 
t E [0, 1]. Denote the sub-curve of C corresponding to [a, b] C [0, 1] as C[a,b]i 
and similarly use C^^.b] for Ci. Denote total curvature as a function TkC )- 



3.1. Total curvatures of piecewise curves 



Definition 3.1 (Exterior angles of PL curves). [IB] The exterior angle be- 
tween two or iented lin e seg ments Pm _iP,„ and Pr„P,„^i , is th e angle between the 
extension of Pm-iPm and PmPm+i, as shown in Figure 2(a) Let the measure of 
the exterior angle to be a 



satisfying: 

< a„ 



< TT. 



This definition naturally generalizes to any two vectors, vi and V2, by joining these 
vectors at their initial points, while denoting the measure between them as ri{vi, V2), 
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as indicated in Figure 
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(a) Orientated lines (b) Vectors 

Fig. 2. An exterior angle 



The concept of exterior angle is used to unify the concept of total curvature for 
curves that are PL or differentiable. 

Definition 3.2 (Total curvatures of PL curves). [TBI The total curvature of 
a PL curve, is the sum of the exterior angles. 

Definition 3.3 (Total curvatures of curves). [1^ The curvature of a 
curve C(t) parametrized on [a, 6] is given by 

|C"(t)xC"(0|| 



Kit) 



t e [a, 6] 



(3.1) 



WC'itW 

Its total curvature is the integral: 1^(^)1 dt. 

Definition 3.4 (Exterior angles of piecevifise curves). For a piecewise 
curve "f{t), define the exterior angle at some ti to be the exterior angle formed by 
7'(ti— ) and y{ti+) where 

liU) - liU - h) 



I'iU-) = lim 



and 



i{U+) = lim 



7(i, + h)^ -fiU) 



h^O h 

Definition 3.5. [Total curvatures of piecewise curves] Suppose that a piece- 
wise curve 4>{t) (regular at the points) is not at finitely many parameters 
ti, - ■ ■ ,tn- Denote the sum of the total curvatures of all the sub-curves as T^i, 
and the sum of exterior angles at ti , • • • , t„ as Tk,2 ■ Then the total curvature of 0(t) 
is T„i -f T„2- 



''This is similarly defined in a recent paper [?]■ 
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3.2. Definitions of convergence 

Definition 3.6. We say that {C^}^ converges to C in parametric measure distance 
if for any e > 0, there exists an integer N such that maxtg[o,i] \Ci{t) — C{t)\ < e for 
all i>N. 

Remark 3.7. For compact curves, this convergence in parametric measure distance 
is equivalent to pointwise convergence. 

Definition 3.8. [ID] Let X and Y be two non-empty subsets of a metric space 
{M,d). We define their Hausdorff distance i.i{X,Y) by 

maxjsup inf d{x,y)^sup inf d{x,y)}. 

Remark 3.9. By the definition of Hausdorff distance, the pointwise convergence 
implies the convergence in Hausdorff distance. 

Definition 3.10. We say that {Ci}J° converges to C in total curvature if for any 
e > 0, there exits an integer N such that |r„(C,) - T«(C)| < e for aU i > iV. We 
designate this property as convergence in total curvature. 

Definition 3.11. We say that {Ci}f^ uniformly converges to C in total curvature 
if for any [^1,^2] C [0, 1] and Ve > 0, there exits an integer N such that whenever 
i > N, |r„{Ci[t^_t2]) ~ '^«;(C[tj,t2])l < £■ We designate this property as uniform 
convergence in total curvature. 

Remark 3.12. Uniform convergence in total curvature implies convergence in total 
curvature. But the converse is not true. 

4. Isotopic Convergence of Inscribed PL Curves 

We will use the concept of PL inscribed curves as previously defined |16j . 

Definition 4.1. A closed PL curve L with vertices vi,V2,--- ,Vm is said to be 
inscribed in curve C{t) if there is a sequence {ij}™ of parameter values such that 
Vi = C{tj) for j — 1,2, ■ ■ ■ ,m. We parametrize L over [0, 1], denoted as L{t), by 

H^j) = Vj for j = 0,1, - ■ ■ ,m 

and L(t) interpolates linearly between vertices. 

The previously established results [Ml Theorem 2.2] and [24l Proposition 3.1] 
showed that a sequence of finer and finer inscribed PL curves will converge in total 
curvature. The uniform convergence in total curvature follows easily. For the sake 
of completeness, we present the proof here. 

Lemma 4.2. For a piecewise curve j{t) parametrized on [0, 1] (which is regular 
at all points), a sequence {Li}'^ of inscribed PL curves can he chosen such that 
{Li}'^ pointwise converges to 7 and uniformly converges to 7 in total curvature. 
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Proof. We first take the end points 7(to) = 7(0) and 7(t„) = 7(1)- And 
then selecl|^ the points where 7 fails to be . Denoted these points as 
{7(^0), 7(^1), •■ ■ : 7(in-i)j 7(^n)}- We then compute midpoints: ) for j G 

{0, 1, . . . , n — 1} to form L2 which is determined by vertices: 

{7(M,7(^),7(ti),---,7(<n-i),7(^^^^^),7(in)}. 

Continuing this process, we obtain a sequence {L^}^ of inscribed PL curves. 

Suppose the set of vertices of Li is {wi.fc = 7(^1. fe)}, for some finitely many pa- 
rameter values ii^fe. Use uniform parametrization [I9] for Li such that Vi^k = Li(ti^k), 
and points between each pair of consecutive vertices are interpolated linearly. Note 
first that this process implies that {Li}f^ pointwise converges to C. For the uniform 
convergence in total curvature, consider the following: 

(1) Consider each tj where 7 fails to be C^. Denote the parameters of two vertices 
of Li adjacent to Li{tj) as tj^ and t*2. Note that limi_j.oo t^ji — hnii^co — ^j- 
This implies that the slope of Li(t'j^)Li{tj) and the slope of Li{t'j2)Li{tj) go to 
j'{tj—) and 7'(tj+) respectively. This shows that 



lim i7(L,(i;.i)L,(t,),L,(<}2)i.fe)) =r;(7'fe-),7'(i,+))- 

(2) Consider a pair of parameters ti i and 2 of any two consecutive vertices of Li. 
Denote the corresponding PL curve as Li[f. ^ ^j. Note that the corresponding 
sub-curve 7i[(. ^^j. of 7 is C^, since the parameters where 7 is not have been 
selected as vertices. Provided that hmi_>oo — ^2,1! = 0, the proof of Theorem 
[ini Theorem 2.2] shows that T^{Li[f. ^ f. — > T^{j[t- ^.t - 2]) as i — 00. 



By Definition 3.5 the above (1) and (2) together imply that the total curvature 
of any sub-curve of Li converges to the corresponding sub-curve of 7 (either 
or just piecewise C^), which is the uniform convergence in total curvature (Defini- 



tion 3.11). □ 



Since uniform convergence in total curvature implies convergence in total cur- 



vature (Definition 3.11), the corollary below follows immediately. 



Corollary 4.3. [16, Theorem 2.2] [24l Proposition 3.1] For C, a sequence {ii}i° 
of inscribed PL curves can he chosen such that {Li^f converges to C pointwise and 
in total curvature. 

Theorem 4.4 (Fenchel's Theorem). [TB] The total curvature of a closed curve 
is at least 27r, with equality holding if and only if the curve is convex. 



''Acute readers may find later that this choice of points is sufficient for this lemma, but not 
necessary. This choice is for ease of exposition. 
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Lemma 4.5. Denote the plane normal to C at some to £ (0, 1) as Il(to). Consider 
two sub-curves C[to-u] o,i^d C[to+v] for some u £ (0,to) o.'f^d v G (io,l)- If both 
Tf,{C[to-u]) < f and Tf,{C[tg+y]) < then these two sub-curves C[to-u] and C[t^+^] 
are separated by n(to) except at C(to)- 



Proof. Denote the point C{to) as a. Suppose that the conclusion is false, then either 
C[to-u] or C[to-i-v] intersects Il(to) other than at a. Assume without loss of generality 
that C[ta-i-v] nn(to) contains another point, denoted as b. Then the sub-curve C[to+v] 
and the line segment ab form a closed curve C[ta-i-v] U ab. So Tf^(C[to+v] U ab) > 2i: 
by Theorem |4.4[ 

Denote the exterior angles at a and 6 as a and jS respectively. Then a = f since 



TT 

2 



n(to) is normal to C'(to)- By Definition |3.1[ /? < tt. By Definition 3.5 we have 

TK{C[t„+v] U ab) = T,,{C[to+v]) + a + ^ < Ti^{C[ta+v] 

So 



Therefore 



> 



which is a contradiction. 



□ 



Theorem |4.6| below is restricted to "inscribed PL curves" . The general theo- 
rem of "piecewise curves, either inscribed or not" will be established later in 
Theorem ES 

Theorem 4.6. For any sequence {Li\f of inscribed PL curves that pointwise con- 
verges to C and uniformly converges to C in total curvature, a positive integer N 
can be found as below such that for all i > N , Li is ambient isotopic to C. 

Proof. For C, there is a non-singular tubular surfac^of radius r |14'. 

Pointwise convergence and the uniform convergence in total curvature imply 
that there exists a positive integer N such that for an arbitrary i > N: 

(1) The PL curve Li lies inside of the tubular surface of radius r; and 

(2) Denote the set of vertices of Li as {wj}"=o- Suppose the sub-curve of C between 
two arbitrary consecutive vertices vj and vj+i as Aj, for j = 0, . . . ,n — 1. Then 
since the total curvature of vjVj^l is 0, the total curvature of Aj can be less 
than |. 



'^We use the terminology of tubular surface as generalization from the recent usage 1141 regarding 
the classically defined pipe surface [17| . 
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Lemma 4.5 implies that all such sub-curves Aj are separated by normal planes 
except the connection points. The facts about fitting inside a tubular surface and 
separation by normal planes provide a sufficient condition |14j for Li being ambient 
isotopic to C. □ 



Remark 4.7. The paper [T3] provides the computation of the radius r only for 
rational spline curves. However, the method of computing r is similar for other 
compact, regular, C^, and simple curves, that is, setting 

r < mm{ , —T—,rend}, 

^max ^ 

where Kmax is the maximum of the curvatures, d„ij:n is the minimum separation 
distance, and rend is the maximal radius around the end points that does not yield 
self-intersections . 



5. Pointwise Convergence 

Pointwise convergence provides a lower bound of the total curvatures of approx- 



imants (Theorem 5.3). The proof relies upon showing this for PL curves first 
(Lemma |5.2[ ). The technique used here is the well known "2D push" 3 . It is suf- 
ficient here to consider a specialized type of push, designated, below, as a median 
push. 

Definition 5.1. Assume that triangle AABC has non-collinear vertices A,B and 
C. Push a vertex, say B, along the corresponding median of the triangle to the 
midpoint of the side AC. We call this specific kind of "2D push" , a median push. 

Lemma 5.2. Let {Li}°^^ be a sequence of PL curves parametrized on [0, 1] and L 
be a PL curve parametrized on [0, 1]. If {Li}°^^ pointwise converges to L, then for 
Ve > 0, there exists an integer N such that Ti^{Li) > Tf^{L) — e for all i > N. 

Proof. For an arbitrary vertex v of L, suppose v = L{tv) for some ty e [0, 1]. Let 
By be a closed ball centered at v. Since L is a compact PL curve, we can choose 
the radius of By small enough such that: 

(1) the ball By contains only the single vertex v of L; and 

(2) it intersects only the two line segments of L which are connected at v. Denote 
these intersections as u = L(tu) and w = L{ty,) for some tmty, G [0, 1]. Then 
u, w and v together form a triangle Auvw. 

Let Ui = Li{tu), Vi — Li{ty) and Wi = Li{tw). Denote the exterior angle of 
the triangle Auvw at v as ri{v), and correspondingly the exterior angle of AuiViWi 
at Vi as ri{vi). (Note that 77(11) is not necessarily equal to the exterior angle of L 
at v. Similarly for rjivi).) By the pointwise convergence we have that the triangle 



ISOTOPIC CONVERGENCE THEOREM 9 



AuiViWi converges to Auvw. So r]{vi) converges to That is, for Ve' > there 
exits an N such that ry(wi) > ri{v) — e' for aU i > N. 

Consider the PL sub-curve of Li lying in By and denote its total curvature as 
T^^{LiDBy). This PL sub-curve of Li can be reduced by median pushes to AuiViWi. 
[TOl Lemma 1.1, Corollary 1.2] implies that T^{Li n By) > 'q{vi). So for i > N, 

T^{L,nBy)>Tj{v)-e'. (5.1) 

Denote the set of vertices of L as V. Then T^{L) = J2vev ''li''^)- Note that 
T^{L^) > E„ev7^«(^» n By). So Inequality [5l] implies that 

> ^«(^* nBy)>Y, V{v) - e'n = T^{L) - e'n 

where n is the number of vertices of L. Let e' = ^, then we complete the prove. □ 

Theorem 5.3. // {Ci}f^ pointwise converges to C, then for Ve > 0, there exits an 
integer N such that Ti^{Ci) > Tf^{C) — e for all i > N. 



Proof. By Lemma 4.2 we can use inscribed PL curves to approximate {Ci}f and 



C, such that the approximations converge pointwise and in total curvature. Then 



apply the Lemma 5.2 to these inscribed PL curves. Since these inscribed PL curves 
converge pointwise and in total curvature to {Ci}f° and C respectively, the desired 
conclusion follows. □ 



6. Uniform convergence in total curvature 

Convergence in total curvature is weaker than uniform convergence in total curva- 
ture. But pointwise convergence and convergence in total curvature together imply 
uniform convergence in total curvature, which is shown by Lemma |6 . 1 1 below . 

Lemma 6.1. If {Ci}f^ converges to a curve C pointwise and in total curvature, 
then {Ci}J° uniformly converges to C in total curvature. 

Proof. Assume not, then there exist a subset [^1,^2] C [0,1] and a t > such 
that for any integer N, there is a i > iV such that \T^{Ci[tj^ ^12]) ~ ^^(^[ti.ta])! > 
that is Tf,{Ci[t^^t2]) > ^K(C[ti,t2]) + t or Tf,{C^[t^^t2]) < ^«;(C[ti,t2]) ~ The latter is 
precluded by Theorem |5.3| Therefore 

T,iC,lt^,t2])>T^iC[t,^t2]) + r. (6.1) 

Consider the sequence of the sub-curves of {Ci}"^ restricted to the complement 
[^1,^2]^ of [ti,t2], and denote it as {Ci[i^,(2]'=}i°- By theorem 
an integer, say M such that for all i > M, 

r,(a,[4,,,,]c) >T„(C[,,,,,]c)-^. (6.2) 



5.3 for 5, there exits 
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Note that r«,(C,) > T^{C,[t^M]) + Tn{C,[t^MY)- So Equations [6l] and ^ imply 
that there is a i > M so that 

T 

r„(Ci) > T^{Ci\t^^t2\) +^fc(C'i[tl,t2]0 > ^«(<^[tl,t2]) + ^K(<^[tl,t2]0 + 2" 

Since C is C^, T,(C[t,,t,]) + T,(C[t,,t,]c) = T,(C). Therefore we get 

which contradicts the convergence in total curvature. □ 



7. Isotopic Convergence 

For a compact curve C, we shall, without loss of generality (Theorem 4.61, 



consider a sequence {Li}'f of PL curves (instead of piecewise curves) as its 
approximation. We shall divide C into finitely many sub-curves, and reduce the 
corresponding sub-curves of Li to line segments, by median pushes, so as to preserve 
isotopic equivalence. The line segments generated by the pushes form a polyline. 
We shall then prove the polyline is ambient isotopic to C. 

To get to the major theorem, we need to first establish some preliminary topo- 
logical results. We use CH{-) to denote the convex hull of a set. 

Lemma 7.1. Let X and Y be compact subspaces of an Euclidean space M**. // 
X n Y = 0, then Y can be subdivided into finitely many subsets, denoted as 
Yi, . . .Yi, . . . Yj-n for some m > 0, such that CH{Yi) n X = for each i. 

Proof. Since X is compact, for Vy G Y, inixex H^; — y|| > 0, and hence 3 an open 
ball By C M.'^ oi y such that By O X ~ 0. Since Y is compact, among these open 
balls, there are finitely many, denoted by By-^ , • • • , By^ such that Y C U"=i ^Vi ■ 
Let Yi = Y n By. for each i — 1, . . . ,m so that 

CH{Yi) = CH{Y nBy^)^ CH{ByJ = By^. 

Thus, for each i, we have CH{Yi) n X = 0. □ 

As we mentioned before, for a simple curve C, there is a non-singular tubular 



surface of radius r (Remark 4.7 1. This surface determines a tubular neighborhood 
of C, denoted as Tq- Denote a sub-curve of C as C'', and the corresponding tubular 
neighborhood of C'' as T'^. 

Lemma 7.2. The compact curve C can be divided into finitely many sub-curves, 
denoted as , . . . ,C'^ , . . . ,C"' for some n > 0, such that 

• T^iC^) < f ,■ and 

• CH{C'') C 
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Proof. By Lemma |7.1[ C can be partitioned into finitely many non-empty sub- 
curves, each which is disjoint from Sr{C). Since C, is also of finite total curvature, 
we can denote these sub-curves as C^, . . . ,C''', . . . ,C" for some n > 0, such that for 
each fc = 1, • • • , n, T^{C^) < f and CH{C'') n Sr{C) = 0. 

Consider C'^ for an arbitrary k — 1, . . . ,n and denote the distinct normal planes 
at the endpoints of C'^ by Hi, 112, respectively. Denote the closed convex subspace of 

that contains C'^ and is bounded by Hi and 112 as H'^. It is clear that CH[C^) C 
H^, but since CH{C^) n Sr{C) = 0, we have that CH{C^) C . □ 

For fc = 1, . . . , n, let \tk-i, tk\ be the subinterval whose image is C^, with corre- 
sponding r*^. Let e be real valued such that 

< e < mm . 

fce{o,...,n} 2 

We extencQ [ife-i,^fe] to [tk^i — e,tk + e], and denote the tubular neighborhood 
corresponding to the extended subinterval as F^, then F^ only intersects F^'+^ and 
F^~^ for each k. 

For a sequence of PL curves {L^}^ converging to C pointwise and in total 
curvature, denote the sub-curve of Li corresponding (with the same parameters) to 

as L^. Denote the end points of by and u^, with the corresponding end 
points of by v^~^ and . 

Lemma 7.3. A large positive integer N can be found such that whenever i > N , 
for each k, we have 

(1) T^{L^) < f ; 

(2) CH{L'^) C F^; and 



(3) |uf — v^\ < ^ and ^(u^ ^u^,C^) < |, where /i(-) refers to the Hausdorff dis- 
tance. 

Proof. Fhc first condition follows from the uniform convergence in total curvature 



(Lemma 6.1 1, and the second and third follow from pointwise convergence. □ 



Now we are ready to reduce each to the segment v!l~^u^ by median pushes. 
In order to prove there is no self-intersection of Li during the pushes, we present 
two lemmas below. The following lemma was established by a recent preprint |llj . 
For the sake of completeness, we give the sketch of the proof here. 

Lemma 7.4 (Non-self-intersection criteria). [TTj Let P = [Pq, Pi, - ■ ■ , Pn) be 

an open PL curve in M.^ . If Tf^{P) — X]j=i '^j then P is simple. 

Proof. Assume to the contrary that P is self-intersecting. Then there must exist 
at least one closed loop. Consider one closed loop. By Fcnchcl's theorem, the total 



"^If C is open and t^—i = or tfe = 1, consider [0, tfc -I- e] or [tk~i — e, 1]. 



12 J. Li, T. J. Peters 



curvature of the closed loop is at least 27r. The total curvature is the sum of the 
exterior angles, among which at most one angle is not counted as an exterior angle 
of P. But an exterior angle is less than tt. So the total curvature of P is at lest 
27r — TT = TT, which is a contradiction. □ 

Milnor [TB] showed the total curvature remains the same or decreases "after" 
deforming a triangle to a line segment, and this can be trivially extended to show 
that the total curvature remains the same or decreases "during" the whole process 



of deforming a triangle to a line segment, as expressed in Lemma 7.5 



Lemma 7.5. If a vertex of a PL curve in undergoes a median push, then the 
total curvatures of new open PL curves formed during the push remain the same or 
decreas^ 

Lemma 7.6. For each k — 1, . . . ,n, use median pushes to reduce L^ to the line seg- 
ment vll^^u^ . Then during these pushes, Li remains simple, and hence the resultant 
PL curve ljfe=i ambient isotopic to the original PL curve Li. 



Proof. Note that the condition ( 1) in Le mma 7.3 implies that T^{L\ ^ \J L'I) < it 
and Ti^{L'^ U L^~^^) < tt. Lemma 7.4 and 7.5 show that the pushed L'^ does not 



intersect its neighbors Lf+i or Lj^"Since CH{L'l) C F^' (Condition (2)), and 
does not intersect for j 7^ fc — 1 or fc + 1, the perturbed L'^ stays inside F^ and 
does not intersect Ll for j 7^ fc — 1 or fc + 1. Then the conclusion follows. □ 

For each k — 1, . . . connecting the end points v^~^ and of , we obtain 
the polyline IJ)!^^ ''^^^^v^ ■ 



Lemma 7.7. The polyline {Jl^i v'' ^v'' is ambient isotopic to [J^^i 



Proof. Perturb to u*, and the line segments move linearly from ^u^ to 
v!l~^v^, and from u^u^^^ to Since - < | and fi{u'^~^uf ,C'^) < 



I (Condition (3) in Lemma 7.3), the perturbation stays inside F^ which has a 



radius r. So during the perturbation, and do not intersect any 

line segments of Li, possibly except their consecutive segments. But note that for 
each fc, Ui, v'' G r^nr^'+^. An easy geometric analysis shows that this restricted area 
of the perturbation precludes the possibility for u'l~^,Ui and u^,v!l^^ intersecting 
their consecutive segments. So the perturbation does not cause intersections, and 
hence preserves the ambient isotopy. □ 



"This holds not only for the median push, but also for any push with a trace lying on the interior 
of a triangle indicated in Definition 15.11 
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Below is the general Isotopic Convergence Theorem, while Theorem 4.6 is a 
special version restricted to "inscribed PL curves" . 

Theorem 7.8 (Isotopic Convergence Theorem). // {Ci}^ converges to C 
pointwise and in total curvature, then there exists an integer N such that Ci is 
ambient isotopic to C for all i > N . 



Proof. For each Ci and e > 0, there exits an inscribed PL curve Li of Ci such 
that Li is sufficiently close (bounded by e) to Ci pointwise and in total curvature 



4.6 to consider, without 



by Lemma 4.2 and ambient isotopic to Ci by Theorem 4.6 Since ambient isotopy 
is an equivalence relation [T^, we now rely on Theorem 
loss of generality, a sequence of PL curves {Li}f instead of {Ci}f^ 

Note that r«(C'=) < f and the polyhne IJ^^^^ 
neighborhood (since CH{C^) C Fc). By the proof of Theorem 4.6 we know that 



these are sufhcient conditions for C being ambient isotopic to IJ)!^]^ w'^ ^v^. By 
the equivalence relation of ambient i sotop y, Lemma |7.7| implies that C is ambient 
isotopic to IJ^^;^ 



and Lemma 



7.6 



to L,. 



further implies that C is ambient isotopic 

□ 



8. Some Conceptual Algorithms and Potential Applications 

The Isotopic Convergence Theorem has both theoretical and practical applications. 
Theoretically, it formulates criteria to show the same knot type in knot theory. 
Practically, it provides rigorous theoretical foundations to extend current algorithms 
in computer graphics and visualization to much richer classes of curves than the 
splines already investigated 



The following are the general procedures derived from our Theorem 4.6 and 



Theorem 7.8 For a specific problem, further algorithmic development will depend 
upon characteristics of the class of curves. If the curve is "nice" in the sense that 
the total curvature and the radius of a tubular surface is easy to compute, then it 
is easy to develop an algorithm. Such "nice" curves include a rational cubic spline 
parameterized by arc length, for which the total curvatures can be easily computed, 
and the radius of a tubular surface can be found according to an existing result |14j . 
Otherwise, for some other curves, the computation of the total curvatures and the 
radius of a tubular neighborhood may be difficult and is beyond the scope of the 
details considered here, even while the theorems provide a broad framework within 
which these subtleties can be considered. 



8.1. Using PL knots to represent smooth knots 

Based on Lemma |4.2| and Theorem |4.6[ a procedure can be designed such that 
it takes a smooth knot as input and picks finitely many points on it to form an 
ambient isotopic PL knot. We call this a PL representation of the smooth knot. 
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Recall that two criteria are sufBcient for the isotopy between a compact 
curve C and its inscribed PL curve C: 

(1) Each sub-curve of C determined by two consecutive vertices of C has a total 
curvature less than ^. 

(2) The PL curve C lies inside of the tubular surface for C with radius r. (This can 
be achieved by making the Hausdorff distance between £ and C less than r.) 

The Outline of Forming PL Representations: 

(1) Select C(0) as the initial vertex of £, denoted as wq. 

(2) Sel[f]e = 0.1. Select|s] ii e [0,1] such that r«(C[o,t,]) = f - £• Let the second 
vertex of £ be C(ii), denoted as vi. 

(3) Similarly pick t2 to obtain V2- Continue until we reach the end point C(l), 
denoted as «„. This process terminates because C is compact. In the end, we 
obtain an £, and sub-curves of C with total curvatures being less than |. 

(4) Verify if the Hausdorff distance between £ and C is less than r; If not, then 
select midpoints: ) for j g {0, 1, . . . , n — 1}, denoted as V 2i+i , to form 
a new inscribed PL curve determined by vertices: 

{wO,Wi,Wl,W|,W2, . . ■ ,Vn-lJV2r^-l,Vn}■ 
{^) Repeat 4 until the Hausdorff distance between £ and C is less than r. This 
process of selecting midpoints implies the pointwise convergence. So this process 
terminates. 

8.2. Testing isotopic convergence 

For a curve C and a sequence {C^}^ of piecewise curves, where {Q}^ 
converges to C pointwise and in total curvature, we shall design a procedure to 
determine a positive integer N such that whenever i > N, Ci is ambient isotopic 
to C. Use Ti^{-) to denote the total curvature, CH{-) the convex hull, and /i(-) the 
Hausdorff distance. 

The Outline of Testing Isotopic Convergence: 

(1) Divide C into sub-curves C*"' for fc = 1, • • • , n such that 

• T^iC") < f ; and 

• CiJ(C'=) C r'=. 

(2) Set i := 1. 

'This e value is not unique. Many others also work. 

sit is not necessary for T'k(C[q ^jj) to be exactly ^ — e. For efficiency, it is fine to end up with a 
value not equal to ^ — e as long as it is less than ^ . This aspect will require a subroutine to be 
developed that will likely vary over the class of curves considered and this detail is beyond the 
scope of the current investigation. 
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(3) Use the above technique to form a PL Representation Li for the piecewise 
curve Ci such that Li is ambient isotopic to Ci. 

(4) Let L^ but the sub-curve of Li corresponding to . Denote the end points 
of L^ as and ttj, and the corresponding end points of as v^~^ and 
i;*^. Verify three criteria: 

• r«(Lf) < f ; 

• CH{L^) C r^; and 

• |uf - w'^l < I and ^{u'y~^u'^,C'') < §. 



(5) If these criteria are satisfied, then let N := i and stop. Otherwise let i := i + 1 
and go to (3). 

We know that {Ci}^ pointwise converges to C and uniformly converges to C in 



total curvature (Lemma 6.1 ), so there exists a finite i such that these three criteria 
are achieved, which means the above process terminates. By Isotopic Convergence 
Theorem, we obtain the ambient isotopy. 



8.3. A potential application for molecular simulations 

It is often of interest to consider geometric models that are perturbed over time. 
For chemical simulations of macro-molecules, the algorithms in high performance 
computing (HPC) environments will produce voluminous numerical data describing 
how the molecule twists and writhes under local chemical and kinetic changes. These 
are reflected in changed co-ordinates of the geometric model, called perturbations. 
To produce a scientifically valid visualization, it is crucial that topological artifacts 
are not introduced by the visual approximations 8 . A primary distinction between 
the approximations created here in Section [8.1| and those based on Taylor's Theorem 
[5] is the expression here of the upper bound for total curvature to be less than f — e. 

An earlier perturbation result is limited to PL curves [2] . Other ambient isotopic 
approximation methods rely upon the curve being a spline |10ll8j . The example in 
the next section needs not be a spline. Our approach provides a more general result 
of sufficient conditions for both approximation and perturbations which do not 
change topological features. Precisely, our Isotopic Convergence Theorem implies 
that as long as the convergence criterion of pointwise convergence and convergence 
in total curvature is satisfied, then ambient isotopy is preserved. 



8.4. A representative example of offset curves 

Offset curves are defined as locus of the points which are at constant distant along 
the normal from the generator curves [13]. It is well-known [22l p. 553 ] that offsets of 
spline curves need not be splines. They are widely used in various applications, and 
the related approximation problems were frequently studied. A literature survey on 
offset curves and surfaces prior to 1992 was conducted by Pham [21], and another 
such survey between 1992 and 1999 was given by Maekawa [13]. Here we show a 
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representative example as a catalyst to ambient isotopic approximations of offset 
curves. 

Let C{t) be a compact, regular, C^, simple, space curve parametrized in [a, 6], 
whose curvature k never equals 1. Then define an offset curve by 

n{t)^C{t)+N{t), 

where N{t) is the normal vector at t, for t G [a, h]. 

For example, let C{t) — (2 cost, 2 sin t) for t £ [0, 27r] be a helix, then it is an 
easy exercise for the reader to verify that the above assumptions of C are satisfied, 
with K = |. Furthermore, it is straightforward to obtain the offset curve Vl{t) = 
(cost, sint, t), which is not a spline. 

We first show that is regular. Let s(t) — \C'{t)\dt be the arc-length of C. 
Then by Frenet-Serrct formulas [S] we have 

n'{t) = c\t) + N'{t) 

ds ^ , „ „^ , ^ ds ^ ds „ 

= -tT + {-kT + tB — = (1 - k)—T + T—B, 
dt ^ ' dt ^ ' dt dt 

where T and B are the unit tangent vector and binormal vector respectively. Since 

T _L B, if (1 - k)^ 9^ then Q'{t) ^ 0. But (1 - k)§ ^ because k ^ 1 and C{t) 

is regular by the assumption. Thus Q{t) is regular. 

Now we define a sequence {ili{t)}°^^ to approximate ^l{t) by setting 

n,it) = cit) + —Nit). 

It is obvious that {rii(t)}^i pointwise converges to Vl{t). For the convergence 
in total curvature, note that \\mi^ao^'i{t) = linii_s-oo (0 = ^"{t), and 

\n'(t)\^Q due to the regularity of Vl{t). Therefore 

^[{t) X vt'ijt) _ n'{t) X n"{t) 

The convergence in total curvature follows. 

Consequently, bjj^the Isotopic Convergence Theorem (Theorem 7.81, we con- 



clude that there exists a positive integer N such that fli{t) is ambient isotopic to 
n{t) whenever i > N. 



9. Conclusion 

We derived the Isotopic Convergence Theorem by topological and geometric tech- 
niques, as motivated by applications for knot theory, computer graphics, visualiza- 
tion and simulations. 

''Since the example satisfies C" and convergence and the paper 1241 shows ambient isotopy 
under C" and convergence, the ambient isotopy for this example also follows from the previous 
result |24| . Here our purpose is to use it as a representation to show how the Isotopic Convergence 
Theorem can be applied. 
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A future research direction could use the Isotopic Convergence Theorem in knot 
classification, since it provides a method to pick finitely many points from a given 
knot, where the set of finitely many points determines the same knot type. 
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